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ENHANCED DENSE GAS FRACTION IN ULTRA-LUMINOUS INFRARED GALAXIES 
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ABSTRACT 

We present a detailed analysis of the relation between infrared luminosity and molecular line lumi- 
nosity, for a variety of molecular transitions, using a sample of 34 nearby galaxies spanning a broad 
range of infrared luminosities (10 10 < Lib, < 10 L@). We show that the power-law index of the 
relation is sensitive to the critical density of the molecular gas tracer used, and that the dominant 
driver in observed molecular line ratios in galaxies is the gas density. As most nearby ultraluminous 
infrared galaxies (ULIRGs) exhibit strong signatures of active galactic nuclei (AGN) in their center, 
we revisit previous claims questioning the reliability of HCN as a probe of the dense gas responsible 
for star formation in the presence of AGN. We find that the enhanced HCN(1-0)/CO(1-0) luminosity 
ratio observed in ULIRGs can be successfully reproduced using numerical models with fixed chemical 
abundances and without AGN- induced chemistry effects. We extend this analysis to a total of ten 
molecular line ratios by combining the following transitions: CO(l — 0), HCO + (l — 0), HCO + (3 — 2), 
HCN(1 — 0), and HCN(3 — 2). Our results suggest that AGNs reside in systems with higher dense 
gas fraction, and that chemistry or other effects associated with their hard radiation field may not 
dominate (NGC 1068 is one exception). Galaxy merger could be the underlying cause of increased 
dense gas fraction and the evolutionary stage of such mergers may be another determinant of the 
HCN/CO luminosity ratio. 

Subject headings: galaxies: evolution — galaxies:fundamcntal parameters and ISM — infrared: galax- 
ies — radio lines: galaxies — submillimeter 
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1. INTRODUCTION 

Over the past 50 years, there have been a number of 
studies relating galaxy star formation rate (SFR) and 
the amount of molecular gas available. iSchmidtJ (|1959D 
proposed a power-law relationship between star forma- 
tion rate volume dens ity and m olecular gas volu me den- 
sity. IKennicuttl (| 19981 ) (also see lKennicuttlll989f) framed 
the problem in terms of surface densities, yielding the 
Kennicutt-Schmidt relation: T,,sfr cx . 

While a power-law index around 1.4 or 1.5 holds 
for the molecular gas traced by CO(l — 0), a power- 
law index of unity was found when re lating SFR to 
the molecular gas traced by H CN(1 — 0) (Solomo n et al.1 
[1991 IGao fc Solomon! l2004bL hereafter GS04b). This 
difference in index is interpreted as a consequence of 
the different critical densities 8 of the molecular trac- 
ers used. Having a lower critical density, CO(l — 0) 
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traces the global molecular gas content, some of which 
may not be involved in the star formation process, 
whereas HCN(1 — 0) traces denser molecular gas (n > 
10 4 cm" 3 ) more cl o sely l inked to star formation . Since 
the ISolomon etafl (|1992f ) and IGao fc Solomonl (|2004aL 
hereafter GS04a) observations, several interpretations re- 
garding the linear I/ir-HCN(1 — 0) relation in galaxies 
have been put forth by both observational arguments and 
numerical models. 

For example, by extending the linear relationship be- 
tween the total infrared luminosity, £ir[= £(8- 1100/im)1, 
and H CN(1 — 0) to Galactic cloud cores, IWu et all 
(2005) added to the original GS04b interpretation by 
framing the relationship in terms of individual star- 
forming dense-gas "units" . As one increases the num- 
ber of dense star-forming clumps in a galaxy, Lir and 
HCN(1 — 0) both increase in lock-step. In this view, the 
only difference between an extreme starburst galaxy and 
Galactic star-forming region is the number of dense star- 
forming units emitting HCN and infrared emission. 

Recent observations may challenge the existence 
of the tight, linear correlation between Lir and 
HCN(l -O) luminosity 9 reported by GS04b. For in- 
stance, IGracia-Carpio et all (puM hereafter GC08) ob- 
tained new HCN and HCO + measurements for 17 lu- 
minous and ultra- luminous infrared galaxies (LIRGs: 
10 11 < L m < 10 12 L Q and ULIRGs: L IR > 10 12 L ). 
They found conflicting HCN(1 — 0) me asurements for 
severa l galaxies that overlap with the [Solomo n et al.l 
(1992) sample. They attribute the differences to cali- 

9 We use the notation L' for molecular line luminosities as they 
are expressed in K km s _1 pc 2 , whereas other luminosities are 
expressed in Lq. 
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bration errors in the older survey, and proceed to com- 
pare their revised Lm— L^cn relation with the original 
one. They find a larger and more significant increase 
of Lin— -^hcn with Lm, meaning that the Lm— LhCN 
power-law would be steeper than an index of unity. 
These authors interpreted the higher Lir/_LJj CN ratio as 
an enhanced star formation efficiency in ULIRGs. 

A large fraction of nearby LIRGs and ULIRGs ap- 
pear to be mergers of gas-rich progenitor galaxies and 
therefore have a high concentration of dense molec- 
ular gas in their center which feeds strong sta rburst 
and AGN activity (see ISanders fc Mirabe]|Tl996l and 
references therein). It has been noted that sources 
with brighter IR lum inosities are both more likely 
to ho s t an AGN (e.g. IVeilleux et all 119951 ; iLutz et all 
119981: IVeilleux et all 119991) . as well as show higher 
HCN(1 — 0)/ CO(1 — 0) molecular line luminosity ra- 
tios (GS04b. IGracia-Carpio et all 120061) . It is tempt- 
ing to associate the increased HCN(1 — 0)/CO(1 — 0) 
luminosity ratio wi t h the presence of AGN. Indeed, 
IGracia-Carpio et all (|2008f ) report evidence for an en- 
hanced abundance of HCN in ULIRGs, and caution that 
care must be taken when converting HCN line luminosi- 
ties into dense gas masses. 

Moreover, some theoretical work suggests that AGN- 
induced processes co uld potentially incre ase the abun- 
dance of HCN (e.g. ILintott fc Viti l2006j) . causing the 
observed Lir-L^ cn relation to have a shallower slope 
than the Lir-L' co relation. In this picture, X-ray emis- 
sion resulting from AGN buried in ULIRGs may en- 
hance the abundance of HCN molecules via an increased 
availability of free electrons that facilitates combina- 
tion with' ions (e.g., HCNH+ + e~ -> HCN + H). 
Some high-resolution HCN observations of nearby galax- 
ies have shown stronger emission of HCN(1 — 0) relative 
to CO(l - 0) in the center of a few LIRGs and ULIRGs 
known to h ost an AGN (e.g. IGracia-Carpio et al.l 120061 : 
lKripsll2007l ). These authors have interpreted this higher 
molecular line ratio in terms of an increased abundance 
of HCN. Others interpret high molecular line ratios 
in terms of the mid-plane pressure, in the sense that 
larger ambient pressure could trap or maintain dense gas 
(|Levine et alj|200l iBlitz fc Rosolowsky||2006D . 

In recent years, new theoretical models suggest a differ- 
ent physical explanation for the observed Lir -molecular 
line luminosity relationships. Specifically, models that 
couple non-local thermodynamic equilibrium (LTE) ra- 
diative transfer calculations wit h models of star-forming 
giant molecular clouds (GMCs) (Krumholz & Thompson 
I2007D and hydro dynamic simulations of di sk galaxies and 
galaxy mergers ([Narayanan et all 12008b?) find that the 
driving force controlling the SFR-L mo i relation in galax- 
ies is the fraction of thermalized molecular gas at a given 
molec ular t ransition. In th e [Krumholz & Thompson 
(|2007l) and iNaravanan et alj ((2008bJ) picture, the SFR 
is controlled by the relation SFRoc n 15 , where n is the 
number density of molecular hydrogen. When the ob- 
served molecular line traces the bulk of the molecular 
gas in the galaxy (e.g. CO(l — 0)), the SFR-L mo i power- 
law index (hereafter "slope" as this relationship is typ- 
ically considered in log-log space) is equivalent to the 
Kennicutt-Schmidt index (i.e. close to 1.5). On the 
other hand, higher critical density molecular line trac- 



ers (e.g. HCN(1 — 0)) trace an increasingly smaller frac- 
tion of the gas in a galaxy, and so the index in the ob- 
served SFR-L mo i pow er-law relation decrease s. These 
theories were tested bv lBussmann et all (|2008f ). who ob- 
served HCN(3 — 2) — a higher critical density tracer 
than HCN(1 - 0) — in the GS04a,b sample of galax- 
ies, and found a SFR-HCN(3 — 2) index less than unity 
(0.72 ±0.08). 

In this paper, we expand on the Bussmann et al. work 
by studying the relationship between five molecular line 
tracers and infrared luminosity. In particular, we em- 
phasize the importance of considering the critical den- 
sity of each molecular line tracer. We present support 
that HCN(1 — 0) is a valid probe of dense molecular 
gas and we revisit the interpretation of the enhanced 
HCN(1 - 0)/CO(1 - 0) luminosity ratio in high IR lu- 
minosity systems. Our results suggest that this feature 
is primarily driven by the increased dense gas fraction in 
galaxies undergoing a merger and/or a strong episode of 
starburst activity; we do not require X-ray driven chem- 
ical abundance effects to explain the observed molecular 
line ratios. Overall, we build a coherent picture where 
the governing parameter is the density distribution of the 
molecular gas. 

We describe our sample of 34 galaxies in fJ5J along with 
the measurements that we use. Section|3]contains our re- 
sults regarding power-law relationships between infrared 
luminosity and various molecular line luminosities (33TTJ) , 
as well as the study of molecular line ratios as a function 
of infrared luminosity ( §3.2|) . The inte rpretation is par- 
tially b ased on the theoretical models oflNaravana n et all 
( 2008b), which we apply to this study in g] We include 
a brief discussion of chemistry effects in [J5] and present 
our conclusions in Throughout this paper we assume 
a flat ACDM cosmological model with (11^,0^, Ho) = 
(0.3, 0.7, 70 km s" 1 Mpc^ 1 ). 

2. METHODS 
2.1. Sample 

Our primary sample consists of 2 9 nearby galax- 
ies w ith HCN(1 — 0) observations (|Gao fc Solomon! 
I2004bl) as well as integrate d optical spectroscopy 
(|Moustakas fc Kennicutl I2006L hereafter MK06). We 
obtained follow-up HCN(3 — 2) observations for 22 of 
these 29 galaxies (10 detections and 12 upper limits) 
usin g the 10 m Heinrich H ertz Submillimeter Telescope 
(see lBussmann et al.ll2008l hereafter B08 for more detail 
on the HCN(3 — 2) observations). In order to expand 
our sample and include additional transitions spanning a 
range of critical densiti es and excitation stat e s, we incor- 
porate the sample from IGracia-Carpio et al. (2008;, here- 
after GC08). These authors include observations of the 
HCN(1 - 0), HCN(3 - 2), HCO+(l - 0), HCO+(3 - 2) 
transitions for 17 galaxies, of which 12 overlap with our 
original sample. Thus, the combined sample comprises 
34 galaxies. We note a distinction between our primary 
sample and that of GC08. While the GC08 sample spans 
the high infrared luminosities (Liu — 10 11,3 — 1O 12 5 L0), 
our sample extends the range down to Lm ~ 10 10 Lq. 

Molecular line and (far-) infrared luminosities are tab- 
ulated in Table [1] When combining our primary sam- 
ple and GC08 sample, we convert all the luminosities 
to common cosmological parameters and luminosity dis- 
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tances. The C0(1 — 0) and HCN(1 — 0) measurements 
were taken from GS04a,b. When available, we substitute 
updated HCN(1 — 0) luminosities from GC08. These au- 
thors published new HCN(1 — 0) luminosities for their 
17 galaxies and found signif icantly differen t value s for 
galaxies overlapping with the ISolomon et al.l (|1992l ) sub- 
sample. These galaxies were also included in the GS04b 
compilation and the new values from GC08 are approx- 
imately a factor of two lower in several cases. GC08 
attributed these discrepancies to observational errors, 
namely the calibration of the receiver used on the 30- 
m IRAM telescope during the earlier observations. 

-^hcn(3-2) vauies were obtained from B08 (see Ap- 
pendix [A| and GC08. Among eight galaxies in com- 
mon, five have consistent values within the uncertainty, 
in which case we use the average (or the maximum up- 
per limit in the case of two upper limits). We adopt 
the GC08 measurements for two non-detections in B08 
(Mrk 231 & IRAS 17208-0014), as well as for NGC 6240, 
which had a discrepant measurement. This galaxy has 
a very broad HCN(3 — 2) emission line, so in this case 
we prefer to use the data from GC08 because of their 
significantly wider bandwidth. 

A few galax i es in our sample were also observed 
by IKrips et al.l (|2008f) . These authors report brighter 
HCN(3 — 2) for five overlapping galaxies, most notably 
for the most IR-luminous galaxies (e.g. Mrk 231). To 
maintain homogeneity in our selected set of observations, 
and to facilitate the comparison with the work of GC08, 
we restrict our analysis to the HCN(3 — 2) luminosities 
from B08 and GC08. However, we anticipate that includ- 
ing the few higher values from IKrips et alj (|2008h would 
either strenghten the trends presented in ij3]or leave them 
unchanged. 

We use the -^Hco+fi-o) an< ^ ^hco+(3-2) measurements 
published in GC08 with a conversion for the luminosity 
distances adopted here. On average, this distance con- 
version changes molecular line luminosities by 0.055 dex 
(min=-0.073 and max=0.16), which is less than the av- 
erage uncertainty on these quantities. 

Observing different molecular transitions with a vari- 
ety of telescopes leads to varying beam sizes. Since most 
sources in the samples we use in this study have not been 
mapped, it is important to consider aperture affects as- 
sociated with the use of varying beam sizes. First, two 
galaxies have a lower limit in HCN(1 — 0) because they 
were not mapped even though they are nearby (NGC 660 
and NGC 2903, see GS04a). For the HCN(3-2) observa- 
tions in B08, the beam size of 30" is sufficient to cover 
the central kpc for galaxies beyond 7 Mpc. We expect 
to detect most of the HCN(3-2) emission in these cases 
since it is the dense nuclear regions of galaxies that are 
responsible for the majority of the emission from molec- 
ular transitions with high critical densities for excitation. 
Meanwhile, targets selected for study in GC08 all lie at 
distances >60 Mpc and therefore we are assured of sam- 
pling the full extent of the high density molecular emis- 
sion. 

2.2. Optical AGN Classification 

The nature of the powering source of galaxies in our 
sample (star formation, AGN, or hybrid i.e. hosting both 
star formation and AGN activity) is determined using 



the o ptical diagnostic diagram known as the BPT dia- 
gram ([Baldwin et al.|[l98lT iVeilleux fc Osterbrocklll987t 
lOsterbrocklll989f) . The optical emission-line ratios used 
in this diagram ([O in] A5007/H/3 and [N n] A6584/Ha) 
probe a combination of the oxygen abundance and ioniza- 
tion parameter of the interstellar medium (ISM) present 
in these galaxies. This provides us with an indirect sig- 
nature of the source of ionizing radiation (young stars 
versus AGN). As shown in Fi gure CD the Sloan Digital 
Sky Survey fSDSS. NrnkeTaXlm M) galaxies ( gray area 
and contours) show a tight excitation seque nce of star- 
formin g galaxies (below the solid line from Kauffmann 
( 2003)), as well as a plume of galaxies that are either hy- 
brid (between the solid and dashe d lines) or AGN (abo ve 
the dashed line adapted from iKewlev et alj (|2001[) ). 
The SDSS emission-line measureme nts are taken from 
the MPA/JHU SPSS team w ebsite 10 (jBrinchmann et all 
12004 iTremonti et all I2004D based on observations of 
galaxies contained in the SDSS data release 7 (DR7; 
lAbazaiian et al.ll2009ft . 

Galaxies from our combined sample are overlaid with 
colored symbols. Among 34 galaxies, 13 are classified as 
star-forming (SF, red circles), 10 galaxies are classified as 
hybrid (SF/AGN, green squares) and 11 galaxies show a 
strong signature of AGN (blue triangles). These numbers 
include five galaxies which are classified but not shown 
on the diagram. Two galaxies (Mrk 231 and NGC 7469) 
have obvious broad emission lines in their optical spec- 
trum indicating a Type 1 (broad-line) AGN, and there- 
fore need not be classified using this narrow-line diag- 
nostic. IRAS 23365+3604 fails the signal-to-noise cut 
for the [N n] A6584 line but occupies the AGN portion 
of the [O i] A6300/Ha diagram. Finally, Arp 299 A & 
Arp 299B are both classified as SF/AGN based on the 
classification of the Arp 299 system. With the exception 
of IRAS 12112+0305 and VII Zw 31, spectra were ob- 
tained with a long-slit drift-scanning technique in order 
to integrate the light spatially (MK06). Spectroscopic 
measurements for IRAS 12112+030 5 and VII Z w 31 were 
obtain ed from IVeilleux et ahl (|1999h and from lWu et alJ 
(1998), respectively. 

We note that dust obscuration can challenge the iden- 
tification of an AGN in the optical spectral range. Suffi- 
cient optically thick material with a large covering frac- 
tion surrounding the nuclear region could potentially 
mask all AGN signatures at these wavelengths. How- 
ever, we identify an AGN in all but one of the most IR- 
bright galaxies (IRAS 17208-0014). The latter is the only 
optically-obscured, potentially MIR-obscured, AGN can- 
didat e sometimes referred to as buried AGN (|Imanishil 
2009). We also note that although this diagnostic in- 
dicates the presence of an AGN, it does not allow to 
quantify its strength and its relative contribution to the 
total infrared light of its host galaxy. For example, the 
well-known ULIRG Arp 220 is thought to be mainly pow- 
ered by star formation (|Lutz et al.lll996D even though i t 



is known to host a Seyfert nucleus (jRieke et alj [l985). 
Furthermore, there has been accumulating evidence that 
most L IRGs and ULIRGs are mainly powered by star- 
bursts (ILutz et al.lll998l:IGenzel et al. (ll998h iFarrah et all 
20031: |yegaeTalT 2008: iNardini et al.ll2008f). except for a 
few notable cases such as Mrk 231 ()Condon et alJll99H 

10 http:/ /www. mpa-garching. mpg.do/SDSS/DR7 
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Fig. 1. — Emission-line diagnostic diagram indicating the posi- 
tion of normal star-forming galaxies (below and to the left of the 
solid curve), AGNs (above and to the right of the dashed curve), 
and galaxies with an admixture of star formation and AGN activ- 
ity (between the solid and dashed curves). The individual galaxies 
in our sample are labeled and plotted using color symbols with 
error bars, with the exception of five galaxies mentioned in the 
text. The different symbols represent the optical spectral type: 
star-forming (SF) (red circles); AGN (blue triangles); and SF/AGN 
(green squares). For reference, the gray-scale and contours (enclos- 
ing 52%, 84%, and 97% of the points) show the locus of emission- 
line galaxies in the SPS S. T he solid and dashed c urves are adapted 
from Kauffmann (2003) and Kcwlcy ct al. (2001), respectively. (A 
color version of this figure is available in the online journal.) 



Spoon et a l. 2007), IRAS 05189-2524 (ISpoon et alJl2007t 
Annus et slT \2007). and NGC 1068 (|Le Floc'h et all 



20011) . These three examples stand out from the rest 



of our sample as they exhibit the largest mid-IR excesses 
(they correspond to the three blue triangles below the 

dotted line in Figure^ 

Among others, (jNardini et al.l [20081 ) report that the 
AGN contribution may dominate in the mid-IR range 
while remaining a small fraction of the total infrared 
light, typically less than 25%. Recent work suggests a 
correlation between the optical spectral classification and 
the AGN contribution to the bolometric luminosity of 
ULIR Gs. Using a sample of 74 ULIRGs, IVeilleux etaTl 
(2009]) find that the AGN contribution ranges from ~ 
15%— 35% among purely star-forming (SF) and LINER 
ULIRGs to ~50 and 75% among Seyfert 2 and Seyfert 

1 ULIRGs, respectively. Individual galaxies may devi- 
ate from this trend as their sample includes a Seyfert 

2 ULIRG with as little as 20% AGN contribution and 
SF-classified ULIRGs with roughly 50% AGN contribu- 
tion. Thus, we stress that although the optical classi- 
fication used here indicates the presence of an AGN, it 
does not allow us to quantify its contribution at infrared 
wavelenghts on a galaxy- by-galaxy basis. Indeed, sev- 
eral systems classified as AGN are also undergoing major 
episodes of star formation. 

2.3. Infrared Luminosity 

We estimate far-infrared [Lfir = L(40 — 500 /im)] 
and total infrared [Ljr = L(8 — 1100 /mi)] luminosities 
of the galaxies in our sa mple using IRAS obs ervations at 
12, 25, 60, and 100 /xm (|Sanders et al.ll2003fh Following 
iMoustakas et al.l (|2006l ). we model the infrared spectral 
energy distribution of each object longward of 100 micron 
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Fig. 2. — FIR/IR luminosity ratio for the MK06 sample of galax- 
ies (gray and black open symbols) and for our combined sample of 
galaxies (filled colored symbols as in Figure [T] also see legend). 
We find an average value of log(L p j r / L j r) = -0.14 (correspond- 
ing to Lju/Lpjn = 1.38, which is close to the value of 1.3 used 
in GC08). Although most of the galaxies in our sample are fairly 
close to that value, there are some noticeable outliers, the most 
striking one being NGC 1068. We stress that our results remain 
unchanged (within the uncertainties) if we adopt FIR instead of 
total IR luminosities in our analysis. (A color version of this figure 
is available in the online journal.) 



using a modified blackbody with dust emissivity propor- 
tional to A^ 1 and dust temperature giv en by the ob- 
served S v (60 nm )/S v (100 fim) flux ratio (| Gordon et al.l 
2000; Bell 2003) . We then integrate numerically over the 
appropriate wavelength range to derive Lfir and Ljr. 

When interpreting galaxy infrared luminosity, one has 
to take into account potential contributions of recent or 
ongoing star formation, older stellar populations, and 
AGN. In the case of LIRGs and ULIRGs, the infrared 
emission from dust heated by powerful starbursts dom- 
inates over that contributed by old stars, so we ignore 
the older stellar population component. It has been sug- 
gested that because the typical AGN IR SED peaks in 
the mid-infrared range (around A = 10 — 20 /im, e.g. 
lElvis et all 119941 ) . the far-infrared portion of the SED 
roughly corresponds to the star- forming component. Wc 
compare the total IR and FIR emission in Figure [5] and 
find a similar ratio (Lm/ Lpm = 1.38) to that reported 
in the literature (1.3, GC08). Most of the galaxies in our 
sample lie within one sigma of the mean value obtained 
for a control sample of galaxies (MK06, small plotting 
symbols in Figure [2]). There are a few outliers, the most 
striking example being the prototype Seyfert 2 galaxy 
NGC 1068, which has the largest MIR-excess and lies 
4.5cr from the mean. Interestingly, even when the frac- 
tional AGN contribution to the mid-IR light dominates 
over the starburst contribution, it still appears to be the 
case that the AGN total IR cont ribution is mostly less 
than 25-30% (jNardini et al.|[200l . 

The larger spread in Ljr/Lfir ratios observed in 
brighter galaxies may be related to a combination of 
varying extinction and AGN contribution. For instance, 
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iVeilleux et al.l (|2009f) suggest an evolutionary sequence 
with three AGN classes for nearby ULIRGs. The first 
class is characterized by small extinctions and large PAH 
equivalent widths. ULIRGs in this class are highly 
starburst-dominated. The second class includes galaxies 
with large extinctions but that are still likely dominated 
by starburst; while galaxies belonging in the third class 
show both small extinctions and PAH equivalent widths 
and have a significant AGN contribution (at least as im- 
portant as the starburst). A detailed analysis is beyond 
the scope of this paper, but we refer the reader to pub- 
lished anal yses utilizing mid-in frared spectroscopy (e.g . 
Lutz et al.lll998t iGenzel et al.lll998t]Spoon et al.l 120071: 
Armus et al.ll2007t iFarrah et al.ll2007fK Here, we simply 
identify the presence of an AGN using optical emission 
lines f tj2.2[) and we do not attempt to quantify its con- 
tribution, given the large uncertainties in doing so. 

Because we do not interpret Lm strictly as a SFR, we 
choose to plot the total infrared luminosity throughout 
this paper. We repeated our analysis using FIR instead 
of IR luminosity. This change has a negligible effect on 
our results (the differences are all well within the uncer- 
tainties, see Tables [5] & El . 

2.4. Numerical Models 

We u tilize the numerical models of iNaravanan et all 
(|2008bfl as a comparison for our observations. We refer 
the reader to this work for details on the models, though 
we summarize the aspects most relevant to the present 
study. 

In an effort to model the obse r ved SF R-L mo i rela- 
tions in galaxies, INaravanan et al.1 (|2008bD coupled 3D 
non-LTE radiative transfer calculations with smoothed 
particle hydrodynamics (SPH) simulations of galaxies in 
evolution. The SPH simulations were calculated uti- 
lizing a modi fied version of the publicly available code 
GADGET-2 (|SpringelH206"5T ). including prescriptions for 
a multi-phase ISM, supernovae pressurization of the ISM, 
and star-formation following a generalized (3D) ver- 
sion o f the Kennicutt-Schmidt law (jSpringel k. Her nquist 
2003). Here, we set the index of the Kennicutt- 
Schmidt relation to 1.5, which has important conse- 
quences in driving the simulated SFR-L mo i relations 
(|Naravanan et alj l2008bl ). We ad ditionally include en - 
ergy feedback from accreting AGN (jSpringel et al.f 2005) , 
though note that it has negligible impact on the simu- 
lated SFR-L mo i relations. 

The molecular line emission properties of the model 
galaxies were extracted using the 3D non-LTE radia - 
tive transfer code, Turtlebeach ((Narayanan et al.ll2006f) . 
Turtlebeach considers both collisional and radiative (de- 
)excitation in determining the level populations of a 
given molecule, and utilizes Monte Carlo methods for 
sampling the spatial and frequency domains. Because 
the hydrodynamic simulations typically have a coarser 
physical resolution than the scale of GMCs, sub-grid 
prescriptions for including GMCs as singular isothermal 
spheres following a Galactic mass spectrum and mass- 
rad ius relation have been impleme nted (for details, please 
see lNaravanan et al.|[200l I2008afl . 

The SPH simulations consist of ^100 galaxies com- 
prised of isolated disks as well as gas- rich, binary, 1:1 
galaxy mergers. The structure o f the g alaxies were ini- 
tialized following the IMo et all (|1998[ ) formalism. In 



order to probe a relatively large dynamic range of 
galaxies, the galaxies were initialized with gas frac- 
tions / g [0.2,0.4,0.8] and total (halo) mass ranging from 
-lxlO 12 M to -4xl0 13 Af Q spaced in 4 mass bins. The 
galaxy mergers were run at a single initial gas fraction 
and mass (/ g =0.4 and M DM = lxlO 12 M ). A key feature 
of these simulations is that they include constant Galac- 
tic abundances for HCN and CO, i.e. no chemistry is 
modeled. These simulations were shown to accurately re- 
cover the observed L m -HCN(1 - 0), L m -CO(l - 0) and 
Lir-CO(3 — 2) relations. Furthermore, they predicted 
a sub-linear Lir-HCN(3 — 2) relation, which was subse- 
quently observed by B08. 

The H CN and CO simulations were taken from the 
study of (Narayan an et al1l2008bf ). while HCO + simula- 
tions of gas-rich galaxy mergers were run specifically for 
this work. 

3. RESULTS 

3.1. Correlation Between Molecular Line and Infrared 
Luminosities 

The molecular gas tracers used in this work span a 
broad range of critical densities (n cr i t = A u i/T u i) varying 
between - 10 33 cm" 3 for CO(l - 0) and ~ 10 7 7 cm" 3 
for HCN(3 - 2) (Figure©. Although the critical density 
of HCN(1 — 0) is an order of magnitude larger than that 
of HCO + (l — 0) it is almost a perfect match to the criti- 
cal density of HCO + (3 — 2). In what follows, we consider 
the values of n cr it at an assumed temperature of 30 K, 
noting that they are roughly constant over a broad range 
of temperatures (20 — 100 K, see Figure©. Even though 
we use n cr it to guide some of our interpretations, we cau- 
tion that molecular gas emissio n can be obs erved in gas 
with densities less than critical (|Evanslll999l ). Substitut- 
ing cri t ical d ensities with effective densities as defined in 
lEvansI (|1999D does not alter our main conclusions. 

Assuming that Lir (or Lfir) provides a good esti- 
mate of SFR and that the molecular line luminosity 
traces the mass of molecular gas above a certain den- 
sity, Lir — L' mol relations (or their surface density equiv- 
alent) are commonly interpreted and/or used to de- 
rive universal SFR prescriptions. In this section, we 
present relationships between galaxy infrared luminosity 
and the luminosity of various molecular lines. We com- 
pute the slope (0) between log(L/fl) and molecular line 
luminosity of CO(l - 0), HCO+(l - 0), HCO+(3 - 2), 
HCN(1 - 0) and HCN(3 - 2) (Figure When avail- 
able, we compute the slope for the combined sample 
(solid lines). However, observations for the HCO + (l — 0) 
and HCO + (3 — 2) transitions are only available for the 
GC08 sample (open symbols). The latter sample com- 
prises fewer galaxies and is biased towards higher in- 
frared luminosities, resulting in derived slopes (dashed 
lines) with much larger uncertainties. Galaxies hosting 
an AGN (blue triangles) occupy the bright end of Lm 
but appear to follow the same relationship as the rest 
of the sample. Molecular line non-detections are shown 
with arrows but are not included when computing the 
slope (f3) or the Pearson correlation coefficient (r). 

All slop es are compu ted using the LINMIX IDL rou- 
tines from iKefivl ([2007) . In these routines, the distribu- 
tions of independent variables are modeled as a mixture 
of Gaussian functions, which allows for greater flexibil- 



G 



Juneau et al. 



^6: 



o ° 



-A A- 



HCN(3-2) 



HC0*(3-2) 
HCN(l-O) 



HCO^l-O) 
C0(3-2) 

CO(l-O) 



10 



100 



T [K] 



Fig. 3. — Critical density, n cr i t = A u i/T u i, for typical molecular 
transitions used to trace molecular gas in galaxies. The critical 
density is fairly insensitive to temperature over a wide range. Note 
that the HCN(1 — 0) and HCO + (3 — 2) transitions have a remark- 
ably close value of critical density, whereas the critical density of 
HCO + (l — 0) is over one order of magnitude smaller than that of 
HCN(1 — 0). Molecul ar data come from the Leiden Atomic and 
Molecular Database (Schoicr ct al. 2005). 



ity when computing the true distributions of these vari- 
ables (i.e. without measurement errors) given the ob- 
servations. Using a Bayesian statistical approach, the 
likelihood distribution function is computed and then 
integrated over the entire data set. This method of- 
fers several advantages compared to other algorithms. 
Namely, it provides likelihood distributions for the val- 
ues of slopes and intercepts while allowing for intrinsic 
scatter (i.e. scatter present in the absence of measure- 
ment errors). Most relevant for this work, this method 
is successful at recovering linear regressions when mea- 
surement errors dominate the scatter and when there is 
a non-negligible number of non-detections in the sam- 
ple. We note that the value of the slope found using 
a simple linear least-squares fitting routine available in 
IDL (LINFIT) is always included within onc-sigma of 
the mean of the distribution of slopes derived by LIN- 
MIX. The mean values of the distribution of slopes are 
shown in Figure [4j We also report the slopes obtained by 
substituting Lm by LpiR (Table 2). The FIR slopes are 
slightly higher on average, but consistent within la. Our 
results agree with previous work in that /?co(i-o) > 1 
(2a) and that Phcn(3-2) < 1 (3c)- The derived slopes 
for HCO+(l - 0), HCN(1 - 0) and HCO+(3 - 2)are sta- 
tistically consistent with each other (and with a slope of 
unity), given the large uncertainties. 

Figure [4] suggests that the proportionality between in- 
frared luminosity log(Lm) and molecular line luminosity 
log(i') varies as a function of critical density. We compile 
the mean and standard deviation of slope distributions 
(j3 ± a) obtained previously and show them as a func- 
tion of the critical density of the corresponding molecular 
transition (Figure [5]) . We find evidence for a shallower 
log(Ljij) — log(L') slope with increasing molecular line 



critical density. We supplement our results with pub- 
lished relations between \og(Lm) and ^og(L' mol ) for the 
following transitions: CO(l - 0) and HCN(l-O) from 
GS04 (open squares), CO(3 — 2) from iNaravanan et al.l 
(2005, filled circle) and, in order of increasing critical 
density, CO(l - 0), C O(2-l), HCO+(l - 0), CS(3 - 2) 
and HCN(1 - 0) f rom IBaan et all ((200S, asterisk sym- 
bols). iBaan et ail (|2008l ) published slopes for the inverse 
relation \log(L' mol ) — log(X/^)]. Thus, the asterisks sym- 
bols are not truly calculated slopes, but rather the inverse 
of the published numbers. These points are included for 
visual reference, but should be considered more uncer- 
tain than the associated error bars. All transitions are 
labeled at the top of the figure. Points corresponding to 
HCN(1 — 0) and HCO + (3-2) are very close given their 
near-identical critical densities. We added small offsets in 
critical density for clarity of the plotting symbols around 
\og(n cr i t ) ~ 3.3 and 6.5. 

Similar results have been found by llono et ail (2009) . 
These authors looked at the reciprocal relation but find 
a slope that varies with critical density in the same 
sense (i.e. the inverse of their slope decreases with n cr i t , 
reaching close to unity for their higher density tracers). 
While they only considered CO(l - 0), CO(3 - 2) and 
HCN(1 — 0), here we extend the trend to a larger num- 
ber of molecular line transitions and, more importantly, 
to a transition with a higher critical density probing the 
regime where the \og(Lui) — log(L') slope is less than 
unity [HCN(3 - 2)]. 

In their coupling of non-LTE radiative transfer calcu- 
lations with hydr o dynam ic galaxy evolution simulations, 
INaravanan et al.l (|2008bD predicted the SFR-L mo i rela- 
tions in local galaxies as a function of increasing molec- 
ular transition which serve as a natural comparison for 
the observations compiled here. We compare the obser- 
vations with the numerical simulations of Narayanan et 
al. in Figure [5] by over-plotting their model predictions 
as the gray shaded regions 11 . The predicted slopes were 
calculated by randomly sampling 25 of the ~100 model 
galaxies in the Narayanan et al. sample, and calculating 
the slope for the subset of galaxies. The shaded region 
is the la dispersion in slopes, and the dashed line is the 
mean. The model SFR-L mo i slope decreases as a function 
of increasing critical density owing to the arguments de- 
scribed in §1. Because higher J-level molecular line tran- 
sitions trace a smaller fraction of the galaxy's molecular 
content, the SFR-L mo i slope decreases from the assumed 
Kennicutt-Schmidt index of 1.5. We find good agree- 
ment between the compiled molecular line observations 
and the theoretical predictions for CO and HCN. 

3.2. Comparison Between High- and Low- Density 
Molecular Tracers 

In this section, we study the possibility that the en- 
hanced HCN(1 — 0)/CO(1 — 0) luminosity ratio actually 
corresponds to an enhanced dense gas fraction. The five 
transitions studied in this work allow us to examine ten 
molecular line luminosity ratios as a function of IR lumi- 
nosity. This set includes published line ratios for which 
the lines involved probe different densities. We intro- 



duce a new line ratio, L'„ 



i-2)/ L I 



HCO+(3-2)/ -^HCN(l-O) 



for which 



11 We note that Narayanan ct al. (2008b) only published the 
model SFR-L mol slopes for CO and HCN. 
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Fig. 4.— Relationship between infrared luminosity and CO(l - 0) (a), HCO+(1-0) (b), HCN(1 - 0) (c), HCO+(3-2) (d), and 
HCN(3 — 2) (e), molecular line luminosities. We show galaxies from our primary sample (filled symbols) as well as from the GC08 sample 
(open symbols). Panels are ordered by increasing molecular line critical density. Our optical AGN classification is shown using the same 
color scheme as in Figure 1 (also see legend). Each panel is labeled with the mean and standard deviation of the corresponding slope (/3) 
distribution and with the linear correlation coefficient r. The gray lines have a slope of unity and are included for visualization purposes. 
(A color version of this figure is available in the online journal.) 



both molecular transitions have nearly the same critical 
density (see Figure [3]). 

When comparing two tracers, we adopt the convention 
of dividing the higher-density (HD) tracer by the lower- 
density (LD) tracer. We quantify the contrast of their 

critical densities as follows: R cr it = ^°s( n cru/ n crit)- Our 
set of line ratios spans 4 orders of magnitude in critical 
density contrast 0.1 < R c nt < 4.4. We quantify the IR 
luminosity dependence of each molecular line ratio by 
a power-law (see Figure [6]). Whenever possible, we fit 
for our combined sample (open and filled symbols, solid 
lines). Otherwise, we fit to the GC08 sample only (open 
symbols, dashed lines). We report the mean and stan- 
dard deviation on the distribution of possible slopes (a) 
in the top of each panel. We also include the critical 
density contrast (R cr it) and the Pearson linear correla- 
tion coefficient (r) in each case. 

Our results are summarized in Table [3] and indicate 
an enhancement in molecular luminosity ratio of high- 



to low-density tracer (L' HD / L' LD ) with increasing IR lu- 
minosity. This trend is observed, with varying degree of 
significance, for all ratios with R cr n > 1. In contrast, 

-^HCO+(3-2)/-^HCN(l-0) — tne 0n ^ ra ti° 01 nnes Wl th 

nearly equal critical densities (R cr it ~ 0) — remains flat 
even at high infrared luminosities. The same trend is 
observed when substituting Ljr by Lfir (see Table [3]). 

The relations are better constrained (> 3cr) for transi- 
tions for which it is possible to use our combined sample 
as it increases the dynamic range in IR luminosity by one 
order of magnitude. Cases where only the GC08 observa- 
tions are available are subject to more uncertainty, and 
would benefit from additional measurements at lower IR 
luminosities to confirm the trend observed here. Nev- 
ertheless, we find a compelling case for interpreting the 
variations in line luminosity ratios in terms of the molec- 
ular gas density distribution. 

In Figure [3 we present a compilation of the values of a 
shown in Figure [5] We identify the points that originate 
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Fig. 5. — Slope of log(Ljn) — log(L') versus molecular line critical densities for the GC08 sample (triangles) and our combined sample 
(stars). Individual slopes are sh own in the previous figu re. We add values from the literature: CO(l — 0) and HCN(1 — 0) from GS04 
(open squares), CO(3 — 2) from Naravanan et al. (2005, filled circle) and, in order of increasing critical density, CO(l — 0), CO(2 — 1), 
HCO + (l — 0), CS(3 — 2) and HCN(1 — 0) from Baan et al. (2008, asterisk symbols). The corresponding molecular line transitions are 
labeled at the top of the figure. Some of the points around log(ra cr it) ~ 3.3 and 6.5 were offset slightly in critical density for clarity of the 
plotting symbols. The shaded regions show the predictions of N08 models for CO (pale gray) and HCN (dark gray) rotational transitions 
from J = 1 — 0toJ = 5 — 4. See text for more detail. (A color version of this figure is available in the online journal.) 



from our combined sample (red stars) and those that cor- 
respond to the GC08 sample (black circles). Although we 
expect other variables such as molecular gas temperature 
and excitation to add to the scatter around each value 
of a, the existence of a significant positive correlation in 
Figure [7] suggests that the molecular gas density distri- 
bution may be the primary physical mechanism driving 
the observed molecular line luminosity ratios in LIRGs 
and ULIRGs. 

Thi s result is consistent with the analysis of llono et alJ 
(2009) who find that the CO(3 — 2) source size is more 
compact than the CO(l — 0) size for ULIRGs suggesting 
that their high FIR (IR) luminosity is linked with them 
having a large amount of dense molecular gas concen- 
trated within their central region. Furthermore, there 
is evidence suggesting that the ISM ambient density 
in ULIRGs is higher by a fa ctor of 100 compare d to 
normal star-forming galaxies (|Solomon et al.lll997f ). A 
number of other studies support the presence of warm 
and dense molecular gas in these IR-b right galaxies (e.g. 
ILahuis et aHl2007t lArmus et alJl2007f >. 

Our findings are also in agr eement with t he Ar p 220 
and NGC 6240 case studies of iGreve et all (|2009t) . who 
report an increased dense gas fraction in these two pro- 
totypical ULIRGs. Using measurements of a large num- 
ber of dense molecular tracers, these authors infer that, 
for these two ULIRGs, most of the molecular gas is in 
a dense phase and that the GMCs mass-size power-law 
is steeper than in normal star-forming galaxies, indicat- 
ing that HCN emission traces denser and more compact 
GMCs compared to HCO+ transitions between the same 



J levels. 

3.3. Possible Causes of Higher Gas Density 

Assuming that variations in temperature and excita- 
tion are smaller than the effect of varying density, a 
positive relation between a and R cr it would imply that 
more IR-luminous galaxies have a larger fraction of dense 
gas than galaxies with smaller IR luminosities (see §4] 
for comparison to radiative transfer models). Several 
mechanisms could create such a situation (which trans- 
lates to a larger average density of the molecular gas). 
For example, superwinds capable of expelling interstel- 
lar mater ial have been observe d in IR-bright starburst 
galaxies (jHeckman et al.1 12000( 1. AGN-d riven outflows 
have been modeled and are predicted to be capable of 
expelling mo lecular gas from galaxies (|Naravanan et akl 
2006. [20081 ). 

Although signatures of atomic gas out flows have been 



repor t ed in a large number of studies (e.g. Hcckman et al 
1990] iVeilleux et al l 119951 : iRupke et al l 120051: lMa~ 
20061: [Spoon fc Holtll2009L for a review, see IVeilleux et al 
(|2005l )). observational evidence for large-scale outflows 
of molecular gas is scarce and very recent. Based on 
their observations of CO(3-2), HCO+(3-2), HCO+(4-3), 
iSakamoto et al.1 ((2009) find P-cygni profiles indicating 
outflows of molecular gas around the two merging nuclei 
of Arp 220. They also report that the HCO + emission is 
more concentrated around the nuclei relative to the CO 
emission, in agreement with the picture presented here, 
where this difference would reflect an enhanced molec- 
ular gas density in these regions. In addition, molecu- 
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Fig. 6. — Molecular line luminosity ratios as a function of infrared luminosity. We compare high- and low-critical density molecular 
lines, \og(L' H D / L' LD ), and label each panel with the best-fit slope (a), the ratio of the critical densities (R cr it = l°g( n ,^it/ n cr?t)) an< ^ 
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The different symbols represent the optical spectral type: star-forming (SF) (red circles); AGN (blue trian- 



gles); and SF/AGN (green squares). Open symbols denote galaxies from the GC08 sample, whereas filled symbols are used for galaxies in 
our primary sample only. (A color version of this figure is available in the online journal.) 
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consistent with being constant with infrared luminosity (a = 0). 
The slope distribution of the points in this figure has a mean and 
standard deviation of 0.10±0.07 (1.5<r away from zero, dotted line). 
(A color version of this figure is available in the online journal.) 



lar gas outflows have recently been identified for the SO 
galaxy NGC 1266 (Alatalo et al., in prep). In this case, 
the outflow may b e associated with t he high-density en- 
shrouded nucleus (|Mouri et alJ ll998). but the powering 
source (AGN or stellar) remains unclear. Future searches 
for molecular outflows in galaxies would be useful to un- 
derstand their energy and gas content balance. 

The gas density distribution could be affected if, for 
example, low-density gas is depleted more easily than 
denser gas during such outflows. In support of this sce- 
nario, lower-density medium is expected to occupy a 
larger volume and would acquire larger outflow veloc- 
ity by conservation of momentum. In addition, the evi- 
dence that ou t flow velocities incre ase with IR luminosity 
(|Martinl[2005t iRupke et al.ll2005f ) is consistent, at least 
qualitatively, with an increase in dense gas fraction in 
more IR-luminous galaxies as they can more easily de- 
plete less dense gas. On the other hand, gas accretion 
and inflows could counteract this effect to a certain ex- 
tent and help replenish the gas content. 

Alternatively, one could posit that galaxies with a 
larger reservoir of dense gas (and/or higher dense gas 
fraction) can sustain higher IR luminosity by allowing 
for more powerful starbursts and AGN activity. Mergers 
of gas-rich galaxies can trigger and fuel both starburst 
and AGN activity and consequently appear as ULIRGs. 
Indeed, most if not all of the ULIRGs in this sample are 
undergoing a galaxy merger. The important change in 
the dynamics of these systems, and in particular, of their 
molecular gas content might be the underlying cause of 
the different gas density distributions and consequently, 



enhanced HD/LD line luminosity ratios. 

Direct comparisons between high- and low-density 
molecular gas tracers to interpret the physical conditions 
of the molec ular ISM in i nfrare d luminous galaxies were 
also used in iBaan et alJ (|2008t ) . These authors propose 
that the gas density distribution reflects the evolutionary 
stage of a nuclear outburst in (U)LIRGs. In their picture, 
the denser gas is depleted more rapidly than the lower 
density gas as the nuclear starburst progresses, and the 
molecular line ratios provide clues about a shift in the 
d ominant heati n g sou rce (UV versus X-ray photons). 

lLoenen et al.l ( 20081 ) mo deled some of the k ey molecu- 
lar line ratios presented in IBaan et alJ (|2008f ) and found 
that mechanical feedback heating is a crucial process to 
explain the low HNC(l-0)/HCN(l-0) ratio observed in 
some systems. They provide additional support for a 
time-dependent model of the physical conditions in the 
nuclear regions of luminous infrared galaxies, according 
to which the molecular ISM switches from a high-density 
(~ 10 5 cm~ 3 ) phase, dominated by stellar heating, to a 
lower-density (~ 10 4 ' 5 cm~ 3 ) phase, dominated by me- 
chanical heating. 

We explore a varying HD /LD line luminosity ratio fur- 
ther in the next section, where we present simulations of 
individual disk galaxies as well as mergers of gas-rich 
galaxies. 

4. MODEL LINE RATIOS 

In order to investigate whether X-ray induced chem- 
istry is necessary to drive the observed trends of line rati o 
with infrared luminosity (e.g. lGracia- Carpio et alJ l2006). 
we compare observed molecular line ratios with fixed- 
abund ance numerical simulations ((Narayana n et all 
I2008H) . 

In the left column of Figure we show the model 
predictions, and compare them directly to the observa- 
tions in the right column. The HCN and CO simulation s 
were taken from the study of (|Naravanan et al.ll2008bl) . 
while HCO + simulations of gas-rich galaxy mergers were 
run specifically for this work. When available, the iso- 
lated disk galaxies include a large dynamic range of gas 
fractions and masses (see § I2.ip . The two galaxy merger 
simulations are identical in all ways except for their feed- 
back implementation. The thick gray curve represents a 
galaxy merger in which 0.5% of the accreted mass energy 
onto the central black hole is re-injected into the sur- 
rounding ISM as thermal energy input, while the merger 
shown by the thin yellow curve does not include AGN 
feedback. Each time step of 5 hr 1 Myr is shown with a 
filled triangle (circle) for the model with(without) black- 
hole feedback. During their peak burst (when the galaxy 
may be most visible as a ULIRG), AGN winds can vary 
both the star for mation history and ga s density profiles in 
the galaxies (e.g. iSpringel et all 2005). Hence, these two 
models effectively serve as two different galaxies prob- 
ing different line ratio-LiR relations. The infrared lumi- 
nosity of the simulated gal axies is es ti mated from their 
known SFR. We use the iKennicuttl (|1998f) conversion 
Lfir = (5.5 x 10 9 ) x SFR and our observed relation 
Ljfj = 1.38 x Lfir found in Figure O 

The molecular line luminosity ratios con- 
sidered here are HCN(1 - 0)/CO(1 - 0), 
HCN(3 - 2)/CO(l - 0), HCO+(3 - 2)/HCO+(l - 0) 



12 



Juneau et al. 



Simulations 



Observations 



0.5 
0.4 
0.3 
0.2 

0.1 
0.0 

0.1 

.05 h 

0.0 
0.6 

0.4 

0.2 
0.0 



▲ merger BH winds 

• merger no winds 

O disk 

ir coalescence 



(a) 



o 
▲ 



o°° a 



O A -JL-* 



HCN(1-0)/C0(l-0) A ▲ AGN (5) 
O • SF 

□ . SF/ AGN 



(c) 



HCN(3-2)/CO(l-0) 



(d) 



AAA 



★ 



o?° 9 



T-r, 1 ^, 



(e) 



HC0+(3-2)/HC0+(l-0) 



A *» . ±#A 




0.8 h 

1 0.6 
j 
\ 

7 0.4 h 

CO 

as 

0.2 - 

0.0 



(g) JlHCN(3-2)/HCN(l-0) 



(h) 



•o 



.0 o 



10.5 11.0 11.5 12.0 12.5 
log (L IR [L ]) 



10.5 11.0 11.5 12.0 12.5 
log (Lb [L e ]) 



Fig. 8. — Several molecular line luminosity ratios as a function of infrared luminosity. The simulation results are shown in the left- 
hand column while the observations for the corresponding ratio are shown on the right. From top to bottom, the luminosity ratios 
are: HCN(1 - 0)/CO(1 - 0)(a,b), HCN(3 - 2)/CO(l - 0)(c,d) , HCO+(3 - 2)/HCO+(l - 0)(e,f), and HCN(3 - 2)/HCN(l - 0)(g,h). This 
figure shows evolutionary tracks of equal-mass merger simulations with (gray line) and without (yellow line) SMBH feedback combined 
with simulations of individual gas-rich disk galaxies (open circles). The tracks start at black-hole coalescence (large star symbol), and each 
time step is marked with a black symbol. Although their dynamical range differs slightly, the models successfully reproduce the spread 
observed in the data. Plotting symbols for the observations are as follows: star-forming galaxies (SF) (red circles); AGN (blue triangles); 
and SF/AGN (green squares). We use open symbols for galaxies that belong to the GC08 sample and filled symbols otherwise. (A color 
version of this figure is available in the online journal.) 



and HCN(3 - 2)/HCN(l - 0). As before, we show the 
galaxies overlapping with the GC08 sample with open 
symbols, whereas filled symbols are for observations 
from our primary data set only. Because galaxies with 
HCO + observations strictly belong to the GC08 sample, 
they occupy the brighter end of the Lir range, which 
is well-sampled by the merger models. The isolated 
gas-rich galaxies (open circles) cover the low-end of the 
infrared luminosity, with the exception of two cases at 
Lm ~ 10 115 — 10 12 i©, which were designed to explore 
extreme conditions that are physically unrealistic at 
low-redshift (e.g. the most massive Mdm ^10 13 Af Q 
galaxies) . 



Generally speaking, there is an excellent correspon- 
dence between the range of the observations and models 
in their molecular line ratios. At large infrared luminosi- 
ties (SFRs), the line ratios tend to increase. This effect 
is simply a manifestation of the fact that the galaxies at 
the high infrared luminosity range are undergoing a star- 
burst event. When the fraction of dense gas in a galaxy 
increases (for example, owing to a merger) so does the 
rate at which stars form. Large amounts of dense gas also 
imply that starbursting galaxies are more easily able to 
excite high critical density tracers (such as various tran- 
sitions of HCN or HCO + ), thus increasing the observed 
HD/LD line ratios. Because the simulations include 



Enhanced Gas Density in ULIRGs 



13 



constant, Galactic-based abundances without chemistry- 
driven variations, the agreement between model and data 
demonstrate that chemistry-driven abundance variations 
are not necessary to produce the observed line ratios. 

We note that generally the dynamic range of the mod- 
eled line ratios seems to match the observations reason- 
ably well. However, we urge caution with a detailed com- 
parison of the models and the data. First, the simula- 
tions were designed to probe a large parameter space 
in gas fractions and galaxy masses. Some galaxies were 
specifically designed to probe relatively extreme condi- 
tions (e.g. initial gas fractions / g =0.8, Mdm ~ 1O 13 M ). 
Consequently, individual model galaxies can be caught 
during a brief snapshot with extreme line ratios, and 
may not exactly map to a particular galaxy from ob- 
served galaxy samples in Figure [8] Second, the galaxies 
in the models were not chosen to precisely mimic the 
relative number of isolated galaxies versus mergers in 
the GS04a,b samples. Thus, the clustering of simulated 
points in Figure [8] may not exactly match those in the 
observations. Third, the Lir in the models was calcu- 
lated using a linear mapping from the SFR. While we 
do not consider AGN contribution to Lib in the simu- 
lations, there is some in the observations. This is likely 
to be the reason why the observations extend to brighter 
Lib. (> 10 12 ' 5 Lq) compared to the simulations, which 
seem to reach a ceiling at Lib ~ 10 122 Lq. 

5. POSSIBLE CHEMISTRY EFFECTS 

It has been suggested that the radiation field associated 
with an AGN influences the abundance of HCN with re- 
spect to other molecules. For example, emitted X-rays 
could cause X-ray dominated regions (XDRs) which have 
different properties than regular photon dominated re- 
gions (PDRs) found around star-forming regions. Pre- 
vious authors claim that conditions ex isting in these 
XDRs affect molecular gas abund ance (|Lintott fc Vitil 
l2006t lAalto et~all[2007[ lKripsll2007D 

On the other hand, feaan et alT(| 200S) used multiple 
molecular line ratios to distinguish between XDR and 
PDR conditions, and found that most of the (U)LIRGs 
in their sample are dominated by PDRs. They used trac- 
ers sensitive to column density (Nh) to distinguish be- 
tween an elevated HCN/CO intensity ratio resulting from 
a high-iVjy PDR or from a \ow-Nn XDR, and found the 
former to be more likely. 

Below we present optical spectral diagnostics associ- 
ated with ISM metallicity and ionization parameter (Fig- 
ure ©. The emission- line ratio [N n] A6584/[0 in] A3727 
has been shown to correlate with the gas-phase oxy- 
gen abundance 12+log(0/i?) while being less affected 
by the pres ence of an AGN than other metallicity 
diagn ostics (jKewlev &; Dopital I2002t IKewlev &; Ellison! 
21 HISS. Meanwhile, the emission-line ratio O32 — 
log([0 ill] A5007/[O III] A3727) traces the ionization pa- 
rameter. The correspondence between the observed ratio 
and the ionization depends slightly on the metallicity so 
the two panels of Figure [9] should be interpreted in con- 
junction. 

Most galaxies are scattered without obvious trends in 
the HCN(1-0)/CO(1-0)-[N 11] A6584/[0 m] A3727 
and HCN(l-0)/CO(l-0)-O 32 planes. The correla- 
tions are very weak or non-existent, especially if we ex- 
clude the two outliers (IRAS 05189-2524 and the more 



extreme NGC 1068). In the first panel, the correlation 
coefficient r = 0.42 drops to r = 0.24 if we exclude 
NGC 1068 and down further to r = 0.15 if we exclude 
both outliers mentioned above. In the second panel, the 
correlation coefficient drops from r = 0.26 to r — —0.09 
if we exclude both outliers. NGC 1068 (and possibly 
IRAS 05189-2524) may be a special case where unusual 
ionization or chemistry effects could play a significant 
role in producing the observed HCN(1 — 0)/CO(1 — 0) 
luminosity ratios. 

Our result is in agreement with findings of lUsero et alJ 
(2004) who conclude that the circumnuclear region of 
NGC 1068 is effectively a giant XDR. These authors 
used a combination of single dish and interferometry 
data of several molecular species to rule out alterna- 
tive ex^hmatiojiB fort he observed high HCN/CO ra- 
tio (|Tacconi et all 119941 ). NGC 1068 is also observed to 
have the largest MIR-excess (it lies 4.5cr from the mean 
on Figure [5]), a regime where MIR-pumping can pro- 
mote HC N(1 — 0) emission by exciting a bending mode 
at 14 jum (jAalto et al.lll995fl . 

While the optical emission line ratios used here may 
give a good indication of the metallicity and ionization 
state of the diffuse ISM, it is not clear how representa- 
tive they are of the conditions present in denser molec- 
ular gas. Therefore, it may not be surprising that we 
do not see obvious trends in either panel. This caveat 
is especially important if young stars (starburst) are the 
main source of ionization because UV photons do not 
penetrate deeply into the dense molecular. On the other 
hand, X-ray photons can penetrate through much larger 
column densities of material, so it is plausible that the 
extreme conditions triggered by an AGN would occur 
over a large enough volume to effect both the diffuse and 
dense gas phases. If this were the case, we would expect a 
larger HCN(1 — 0)/CO(1 — 0) luminosity ratio to couple 
with larger values of O32 because an increased ioniza- 
tion means more free electrons, which can accelerate the 
production of HCN molecules (as discussed in fjl|) . 

Overall, our results indicate that cases with a genuine 
abundance change in HCN may exist but are the excep- 
tion rather than the rule. 

6. CONCLUSIONS 

Using a sample of 34 nearby infrared luminous galaxies 
(10 10 < Lib < 10 12 5 Lq), we characterized the infrared 
luminosity dependence of various molecular gas tracers. 

• The presence of AGN was assessed using the op- 
tical BPT diagram. In agreement with previous 
publications, we find a more frequent occurrence 
of AGN in more IR- luminous galaxies. This may 
be related to the availability of larger amounts of 
dense gas during the mergers of gas-rich galaxies. 

• The molecular transitions used, CO(l — 0). 
HGO+(1-0), HCN(l-O), HCO+(3-2) and 
HCN(3 — 2), span 4 orders of magnitude in critical 
density. We find that the relationship between 
(F)IR luminosity and molecular line luminosity 
L mo i is shallower for transitions with higher n cr n. 
This trend is in agreement with theoretical models 
of N08 and KT07 and can be explained by the 
varying degree of thermalization of the gas giving 
rise to molecular line emission. 
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Fig. 9. — HCN(1 — 0)/CO(1 — 0) luminosity ratio as a function of nebular emission-line diagnostics of the diffuse ISM metallicity (Left) 
and ionization parameter (Right). There is no obvious trend except for two outliers for which the high HCN(1 — 0)/CO(1 — 0) luminosity 
ratio could potentially be associated with enhanced metallicity and/or ionization in the ISM. Plotting symbols represent the optical spectral 
type: star-forming (SF) (red circles); AGN (blue triangles); and SF/AGN (green squares). We use open symbols for galaxies that belong 
to the GC08 sample and filled symbols otherwise. (A color version of this figure is available in the online journal.) 



• Trends of molecular line ratios with Lir are con- 
sistent with an increased molecular gas density 
in more IR-bright galaxies. T his result agrees 
with t he pi cture pres e nted by iGao fc Solomon! 
(2004b]) and IWu et all (p005h . When comparing 
high-density and low-density molecular gas trac- 
ers, we observe an increase in their luminosity ra- 
tio (Lhd/ L>ld) with increasing infrared luminos- 
ity. Interestingly, this trend vanishes when com- 
paring two tracers with nearly equal critical den- 
sity (HCN(1 - 0)and HCO+(3 - 2)). We infer that 
the main driver of the enhanced HD/LD luminos- 
ity ratio is the molecular gas density distribution, 
in agreement with the observations that ULIRGs 
host large reservoirs of dense molecular gas in their 
central regions. 

• We compare our observed molecular line ratios with 
theoretical values obtained from a set of galaxy 
models. We consider SPH simulations of two 
galaxy mergers as well as individual gas-rich galax- 
ies. With constant Galactic abundances, the mod- 
els successfully produce an enhanced HD/LD lumi- 
nosity ratios at brighter infrared luminosity. This 
provides additional support for a higher molecular 
gas density in galaxies that are extremely gas-rich 
or undergoing gas-rich mergers. Indeed, this result 
suggests that AGN-induced chemistry (or other) 
effects may not be necessary to reproduce the ob- 
servations, but rather that AGNs are more likely to 
reside in galaxies that are very gas-rich and/or ex- 
periencing a merger. The simulations also demon- 
strate important variations in molecular luminos- 
ity ratio with the evolutionary stage of the merg- 
ers. These variations may be connected with the 
dynamics of the molecular gas during mergers (in- 
flows, outflows, gas compression, etc.). 



• We investigate possible chemistry effects on the 
well-known HCN(1 - 0)/CO(1 - 0) luminosity ra- 
tio using common optical emission line diagnostics 
of the ISM ionization parameter and metallicity. 
With the exception of one of two outliers, we do not 
observe significant correlations between the molec- 
ular line ratio and these properties. Decoupling be- 
tween the dense molecular gas traced by CO(l — 0) 
and HCN(1 — 0) and the more diffuse gas that 
traces the optical nebular line may be responsi- 
ble for the absence of correlation. However, we 
note that NGC 1068 shows extreme conditions in 
the sense of having the largest optically-measured 
ionization parameter as well as the largest MIR- 
excess. This outlier is thus more subject to chem- 
istry and/or MIR-pumping effects quoted in the lit- 
erature. We emphasize that NGC 1068 has prop- 
erties that are very distinct from the rest of our 
sample, and thus may not be representative of its 
class. 



• More high-density molecular line observations 
would be extremely beneficial to confirm the trends 
outlined in this work. Only half of our combined 
sample of 34 galaxies have data for all five molec- 
ular lines used here. In particular, future observa- 
tions should target galaxies with fainter IR lumi- 
nosities to test whether the results presented here 
extend into the lower luminosity regime. Higher- 
resolution studies of high-density tracers may com- 
plement this work by allowing for a more detailed 
analysis of the density and excitation structure of 
the molecular gas within galaxies. Such studies 
could provide information on the local influence of 
AGN, whereas in this work we probe the global in- 
fluence of AGN. 
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• Overall, our results support that HCN(1 — 0) is a 
valid tracer of dense molecular gas in galaxies even 
in the presence of AGN. We expect scatter in the 
relationship due to variations in temperature and 
possible radiative (de-)excitation of this transition. 
Also, the dense-to-total molecular gas fraction is 
expected to differ from galaxy to galaxy, especially 
in systems undergoing significant mergers. 
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TABLE 1 

Infrared, far-infrared and molecular line luminosities for our combined sample of 34 galaxies. The optical classification is 

included in the last column. 





Source 


Dl 


Lir 


Lfir 


; 

-^CO(l-O) 


j 1 a 

1J HCN(l-0) 


r / 

^HCN(3-2) 


Y i 

HCO+(1-0) 


Y i 

HCO+(3-2) 




Class 




(Mpc) 


log(L Q ) 


l°g(£©) 


log(L') 


log(L') 


log(L') 


log(L') 


log(L') 




Tf~* 1 i^OQ 

1C lozo 


84 i 6 


11.64 


± 0.07 


11.4S ± U.U/ 


1 a 1 A _i_ n m 
1U.14 ± U.IU 


8.95 


it 0.10 








or? 
or 


Mpri con 


37 i 3 


11.04 


± 0.09 


1 n riQ _L n no 
lU.yo ± U.Uy 


n q q _i_ n 1 n 

y.oo ± U.IU 


7.97 


it 0.10 


T 1 c _i_ n no 
(.10 ± U.Uy 






or /ACjIN 


Mnri can 
INLiO OOU 


lo ± 1 


10.50 


± 0.11 


1 n Q7 i n 11 
lU.o/ ± U.ll 


Q on i n in 
o.oU ± U.IU 


> 


7.35 


R(\ I fi 11 

d.du ± u.ii 






or / AOIN 


WPP (iQC. 


1 3^ 4- 7 
loo IC / 


11.62 


it 0.05 


1 1 AS -1- fl O 1 ^ 


q no _i_ n 1 n 
y.yo in u.iu 


8.29 


it 0.09 


/ 7 7t; 


S AA 4- fl OR 
O.^^l It U.UD 


< t .o*± 


or 


IN vjU lUUO 


16 i 1 


11.41 


i 0.10 


ii r\A -U n i n 

_L J. .VJ'-k □! U.1U 


q qi _i_ n 1 n 


8.55 


it 0.10 


7 c;o _|_ n HQ 
i .oo u.uy 






AGN 


NGC 1144 


119 ± 7 


11. oy 


It U.UO 


11.27 ± 0.06 


10.06 =b 0.10 


o.4o 


it U.IU 








SF 


NGC 1614 


67 ± 5 


11.60 


it 0.08 


11.40 ± 0.08 


9.45 ± 0.10 


8.15 


it 0.10 


< 7.36 






SF/AGN 


n^lRQ 9^9 /L b 


1 sn -t- 7 
-Lou in i 


12.08 


± 0.04 


11.04 HI U.U4: 


y.o/ zc u.iu 


8.84 


± 0.10 








AGN 


hpp 914fi 

IN VJV ZrJ_1U 


17 it 2 


11.08 


± 0.10 


i n 4- n i n 

lu.yo m u.iu 


n oa 4- o 1 
y.Zr'i m u.iu 


8.12 


± 0.10 


U.ll ^ u.iu 






SF/ AON 

or / in.vjHN 


Arp 55 


170 ± 7 


11.66 


± 0.04 


11.56 it 0.04 


10.14 ± 0.10 


8.69 


± 0.04 


< 8.17 


8.65 ± 0.06 


8.11 it 0.12 


SF 


NGC 2903 


8 ± 1 


10.23 


it 0.18 


10.08 ± 0.18 


8.68 it 0.10 


> 


7.27 


5.95 it 0.10 






SF 


UGC 05101 


170 ± 7 


11.93 


± 0.04 


11.84 ± 0.04 


9.76 ± 0.10 


8.98 


± 0.04 


< 7.87 


8.70 it 0.07 


8.45 ± 0.04 


AGN 


NGC 3079 


19 ± 1 


10.76 


± 0.09 


10.65 ± 0.09 


9.54 ± 0.10 


8.16 


± 0.10 


7.20 it 0.08 






SF/AGN 


NGC 3628 


11 ± 1 


10.36 


it 0.09 


10.24 ± 0.09 


9.24 ± 0.10 


7.77 


it 0.10 


6.56 it 0.10 






SF 


NGC 3893 


15 ± 1 


10.14 


it 0.09 


10.00 it 0.09 


8.71 ± 0.10 


7.46 


it 0.10 


< 6.00 






SF 


NGC 4414 


17 ± 


10.54 


it 0.02 


10.40 it 0.02 


9.22 ± 0.10 


7.76 


it 0.10 


< 5.70 






SF 


Arp 193 


106 ± 7 


11.64 


± 0.06 


11.55 ± 0.06 


9.72 ± 0.10 


8.49 


it 0.03 


7.88 it 0.04 d 


8.70 ± 0.03 


8.23 ± 0.03 


AGN 


NGC 5194 


8 ± 1 


10.33 


it 0.13 


10.20 it 0.13 


9.13 it 0.10 


7.54 


it 0.10 








SF 


UGC 08696 


165 ± 7 


12.09 


± 0.04 


11.98 ± 0.04 


9.89 it 0.10 


8.87 


it 0.05 


8.56 it 0.09 c 


8.88 ± 0.05 


8.53 ± 0.04 


AGN 


Arp 220 


85 ± 6 


12.16 


it 0.07 


12.08 it 0.07 


10.01 it 0.10 


9.23 


it 0.01 


8.82 ± 0.02 d 


8.89 ± 0.02 


8.23 ± 0.02 


AGN 


NGC 6240 


111 ± 7 


11.84 


± 0.06 


11.69 ± 0.06 


10.01 ± 0.10 


8.89 


± 0.04 


8.78 it 0.04 c 


9.15 ± 0.02 


8.45 ± 0.06 


AGN 


17208-0014 b 


187 ± 7 


12.32 


± 0.04 


12.25 it 0.04 


10.24 ± 0.10 


9.27 


it 0.03 


8.87 ± 0.04 c 


9.10 it 0.04 


8.67 it 0.05 


SF 


NGC 6701 


61 ± 5 


11.07 


it 0.08 


10.94 ± 0.08 


9.59 it 0.10 


8.20 


it 0.10 


7.43 ± 0.13 






SF/AGN 


NGC 7130 


71 ± 6 


11.35 


± 0.07 


11.22 ± 0.07 


9.73 ± 0.10 


8.59 


it 0.10 








SF/AGN 


IC 5179 


50 ± 4 


11.18 


it 0.08 


11.05 it 0.08 


9.49 it 0.10 


8.61 


it 0.10 








SF 


NGC 7469 


69 ± 2 


11.59 


it 0.04 


11.39 it 0.04 


9.59 ± 0.10 


8.43 


it 0.02 


7.87 it 0.05 d 


8.59 it 0.01 


7.79 it 0.09 


AGN 


23365+3604 b 


269 ± 6 


12.12 


it 0.02 


11.97 it 0.02 


9.94 ± 0.10 


8.78 


it 0.08 


< 8.47 f 


8.61 it 0.10 




AGN 


NGC 7771 


60 ± 5 


11.35 


it 0.08 


11.23 it 0.08 


9.96 ± 0.10 


8.63 


it 0.01 


7.85 ± 0.03 d 


8.62 ± 0.03 


7.67 it 0.08 


SF 


Mrk 331 


77 ± 6 


11.42 


± 0.07 


11.28 ± 0.07 


9.75 it 0.10 


8.41 


± 0.02 


7.66 it 0.10° 


8.50 ± 0.02 


7.94 ± 0.07 


SF/AGN 


Mrk 231 c 


188 


12.51 


± 0.07 


12.26 ± 0.07 


10.00 it 0.10 


9.09 


it 0.03 


8.49 ± 0.05 e 


9.07 ± 0.04 


8.50 ± 0.06 


AGN 


12112+0305 b ' c 


335 


12.31 


it 0.07 


12.19 ± 0.07 




9.22 


± 0.09 


< 8.90 e 


8.99 ± 0.11 




AGN 


VII Zw 31 c 


241 


11.96 


± 0.07 


11.83 it 0.07 


10.16 ± 0.10 


8.84 


± 0.07 


< 8.57 


8.90 ± 0.06 


8.44 ± 0.10 


SF 


Arp 299A C 


50 ± 4 


11.67 


± 0.07 


11.48 ± 0.07 




8.21 


± 0.02 


7.45 it 0.11° 


8.48 ± 0.02 


7.93 ± 0.04 


SF/AGN 


Arp 299B C 


50 ± 4 


11.56 


± 0.07 


11.29 ± 0.07 




8.02 


± 0.03 


< 7.19 e 


8.22 ± 0.03 


7.55 it 0.07 


SF/AGN 



a When no formal uncertainties are available for published CO(1-0) or HCN(l-O) line luminosities, we estimate the calibration error to be 0.1 dex (20—25%). 
b Full galaxy name starts with IRAS. 

c These five galaxies are not part of our primary sample. Lir and Lfjr were obtained from GC08 after converting their luminosity distance to our assumed 

cosmology ({Qm,^a, Ho} — {0.3, 0.7, 70 km s _1 Mpc -1 }). Because GC08 do not provide errors on Ljr and Lfir, we use the mean error for the rest of 

our sample (~ 0.07 dcx). 

d Average of B08 and GC08 values. 

e Measurement from GC08. 

f Upper limit in both B08 and GC08 surveys. 
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TABLE 2 

Slopes of log(L ir ) - log(L^ o; ) and log(L fir ) — hOG(L' mol ) for five molecular 
lines. The transitions are ordered in increasing critical density (n crit ). For each 
transition, we report the difference between the two values of slopes 
(Pfir-Pir) and we identify the sample used for the calculations. 



Transition 


\og(n cr it[cm 3 ]) 


PlR 


Pfir 


Pfir-Pir 


Sample 


CO( J = 1-0) 


3.3 


1.30 ± 0.17 


1.34 ± 0.16 


0.04 


combined 


HCO+(J = 1-0) 


5.3 


0.99 ± 0.26 


1.07 ± 0.26 


0.08 


GC08 


HCN(J = 1-0) 


6.5 


1.13 ± 0.11 


1.14 ± 0.10 


0.01 


combined 


HCO+(J = 3-2) 


6.6 


0.81 ± 0.21 


0.90 ± 0.19 


0.09 


GC08 


HCN(J = 3-2) 


7.7 


0.70 ± 0.09 


0.71 ± 0.08 


0.01 


combined 



TABLE 3 

Slopes of log(L' hd /L' ld ) - wg(L ir ) and log(L' hd /L' ld ) - 
log(Lfih) for ten molecular line luminosity ratios. For each line 
ratio, we report the contrast in critical density (rcrit), the mean 
and standard deviations of the slope distributions (aj r & ct fir ), 
and the linear correlation coefficients (rj r & tfir.)- 



Line Ratio 




O-IR 


&FIR. 


rm 


TFIR 


HCOgt 


. 2 /HCNi_ 


0.1 


-0.054 ± 0.24 


-0.068 ± 0.22 


-0.05 


-0.08 


HCNi_ 


-o/HCO+_ () 


1.1 


0.32 ± 0.15 


0.36 ± 0.13 


0.60 


0.66 


HCN 3 - 


-2/HCO+ 2 


1.2 


0.32 ± 0.39 


0.38 ± 0.38 


0.37 


0.42 


HCO+_ 2 /HCO+ „ 


1.3 


0.24 ± 0.17 


0.26 ± 0.15 


0.48 


0.53 


HCN 3 - 


. 2 /HCNi_ 


1.3 


0.41 ± 0.11 


0.41 ± 0.12 


0.70 


0.72 


HCO+ 


. /COi_o 


2.0 


0.24 ± 0.20 


0.22 ± 0.19 


0.40 


0.28 


HCN 3 _ 


- 2 /HCO+_ 


2.4 


0.53 ± 0.31 


0.62 ± 0.30 


0.68 


0.75 


HCNi_ 


-o/COi_o 


3.1 


0.23 ± 0.064 


0.22 ± 0.065 


0.53 


0.49 


HCO+ 


. 2 /COi_ 


3.2 


0.48 ± 0.27 


0.48 ± 0.26 


0.60 


0.60 


HCN 3 _ 


- 2 /COi_o 


4.4 


0.70 ± 0.12 


0.68 ± 0.11 


0.87 


0.87 



18 



Juneau et al. 



APPENDIX 
OBSERVED HCN(3-2) SPECTRA 
The Heinrich Hert z Submillimeter Telescope (HHT) HCN(3 — 2) spectra from the survey presented in 



iBussmann et a l. (2008) are displayed in Figure [TOl 
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Fig. 10. — HCN(3-2) spectra for all galaxies observed in B08. The 23 galaxies overlapping with this sample are marked with an 
asterisk. One galaxy was mapped (NGC 253, plus symbol). These observations were combined with the HCN(3-2) data presented in 
IGracia -Carpio et al. ( 2008) as described in £12.11 
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Fig. 10. — Continued. 
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Fig. 10. — Continued. 



